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Abstract 

We study the isotropisation of the homogeneous but anisotropic 
Bianchi class A models in presence of a minimally coupled and massive 
scalar field with or without a perfect fluid. To this end, we use the 
Hamiltonian formalism of Arnowitt, Deser and Misner(ADM) and the 
dynamical systems analysis methods. Our results allow to define three 
kinds of isotropisation called class 1, 2 and 3. We have specifically 
studied the class 1 and obtained some general constraints on scalar- 
tensor theories which are necessary conditions for isotropisation. The 
asymptotical behaviors of the metric functions and potential in the 
neighborhood of isotropy have also been determined when the isotropic 
state is reached sufficiently quickly. We show that the scalar field 
responsible for isotropisation may be quintessent and that the presence 
of curvature favor a late times acceleration and quintessence. Some 
applications are made with the well known exponential law potential 
by using our theoretical results but also by help of numerical analysis. 
The isotropisation process with a power law potential is also explored. 
We think this work represents a framework able to guide some future 
researches on the isotropisation of homogeneous models in scalar-tensor 
theories and we argue by discussing briefly about some recent results 
we have obtained in presence of a non minimally coupled scalar field 
or several scalar fields. 

Please, review the English for style 



1 Introduction 

Two basic elements used to define a cosmology are a geometry represented 
mathematically by a metric and a material contents described mathemati- 
cally by a Lagrangian. In this work, we will consider the isotropisation of the 
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homogeneous but anisotropic cosmologies of Bianchi in presence of a scalar 
field. Let us begin by explaining what a scalar field is and why it seems 
impossible to circumvent in cosmology. 

By definition, a scalar field is a function which at each point of space and 
time associates a number. A good example is the temperature of a room: in 
each one of its points, one can associate a number defining a temperature. 
Another example is the gravitational potential 4> outside a mass M . These 
fields abound in particles physics and although they were not detected yet, 
one can advance some arguments making their prediction natural. Hence, 
the best experimentally established theories rest on vectors and tensors fields: 

• the electromagnetic interaction is described by a neutral and massless 
vectors field corresponding to a photon without any charge and an 
interaction of infinite range. 

• the weak interaction is described by a massive vectors field representing 
the W and Z bosons. It means that it is with short range and unstable: 
these particles disintegrate in pairs of other particles. 

• the strong interaction is described by a massive vectors field of gluons. 

• the gravitational interaction is described by a tensors field correspond- 
ing to graviton. 

All these vectors and tensors fields correspond to particles with integer spin 
being respectively Ifi^ and 2%, called bosons^. It then seems natural to won- 
der what a boson of spin would represent: the answer is "a scalar field". 
Probably the most known is the field of Higgs-Englert, introduced in order to 
solve a renormalisation problem[T] involved by the weak force: rather than 
to impose a mass on bosons propagating the weak force, it is supposed that 
the vectors field representing them interacts with a massive scalar field (f), 
i.e. having a potential which describes the interaction with the vectors field. 
The form of this potential would be U{4>) = /c((/>^ — (l)^). This field, which 
one hopes for detection by the LHC in 2008, gives their mass to the particles 
which interact with it. 

Reasons different from those of Higgs can be evoked concerning the presence 
of scalar fields in particles physics. Thus, the supersymmetric theories pos- 
tulate the equality between the fermionic and bosonic degrees of freedom: 
to each boson(of which that of Higgs), a fermion corresponds and vice versa. 
However, this can be done only by adding some degrees of freedom which 
are introduced with help of scalar fields whose potentials U {(p) can be com- 
pletely different from the one above. Generally these scalar fields are called 

is the constant of Planck divided by 2tt 
^Particles with half spin belongs to the family of the fermions. 
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dilaton. The supersymmetric particles could be essential components of the 
dark matter. In particular the lightest of the neutralinos, a state resulting 
from a mixture of several supersymmetric particles, could be the lightest of 
the supersymmetric particles and a candidate for cold dark matter. For a 
popular introduction to the supersymmetry theory, one will refer to the book 
of Gordon Kane(21, "Supersymmetry, unveiling the ultimate laws of Nature". 
If these fields are present in particles physics, it is likely that they populate 
the Universe and thus must, like any material contents, have some cosmolog- 
ical consequences. The first relativistic scalar-tensor theory was the Jordan, 
Brans and Dicke theory^ in the sixties which was issued^ of the ideas 
of Kaluza-Klein0, Diracjl] and Mach. More modern reasons to consider 
scalar fields in cosmology appear at the beginning of the eighties: the ideas 
of Guth on inflation gave birth to scalar fields called inflatons. At that time 
cosmology suffers from many conceptual problems: why the Universe seems 
so flat? How causally disconnected areas at the beginning of times can be 
so similar today? Guth|7| noticed that they would be partially solved if, at 
early times, the Universe undergoes a period of inflation with an exponential 
expansion. For that, the first idea which comes to mind is to introduce a 
cosmological constant. However the observations show that its current value 
would be lO^^'^ time smaller than that predicted at early times: it is the 
cosmological constant problemjHj. A way to solve it is to consider a massive 
scalar field, i.e. having a potential representing the coupling of the field with 
itself and able to play the role of a variable cosmological constant. Since the 
end of the nineties, the presence of this potential found new reasons to exist 
with detection by two independent teams^ JUj of the Universe expansion 
acceleration. One of the most widespread explanations of this phenomenon 
would be the presence of a quintessent scalar field whose density and pres- 
sure are bound by an equation of state similar to that of a perfect fluid and 
whose barotropic index would be negative. It then results a negative scalar 
field pressure which would be at the origin of this new and recent period of 
expansion acceleration. 

With regard to the geometrical description of the Universe, we will be inter- 
ested by the homogeneous but anisotropic cosmological models of Bianchi 
for which the expansion rate depends on the observation direction: it is thus 
an anisotropic Universe on the contrary traditional models of Friedmann- 
Lemaitre- Robertson- Walker (FLRW) where the expansion is the same what- 
ever the observation direction. Let us start with a short history of these 
models. Luigi Bianchi^ was born in Parma on June 18, 1856. It was the 
student of Ulisse Dini and Enrico Batti at the Ecole Normale of Pisa and 
became professor at the University of Pisa in 1886 then the Ecole Normale 

^ These information come from |http : / /www34 .homepage .villanova. edu/ robert . j antzen / 1 
bianchi /bianchi. html^j^papers. 
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director in 1918 until its death in 1928. Its mathematical contributions were 
published in 11 volumes by the Italian Mathematical Union and cover a great 
number of fields. With regard to the riemanian geometry, he is especially 
known for its discovery of the Bianchi identities whose full demonstration 
was given in ^^(it had discovered them first once in an article of 1888 but 
had neglected their importance by giving them in footnote). In 1897, by 
using the results of Lipshitz[T2I and Killing[TB| as well as the Lie^l] the- 
ory of continuous groups, he gave a complete classification of the Riemann 
3-varieties isometrics classes, identified by the roman letters / to IX. At 
the time, neither general relativity, nor special relativity existed yet^. In 
1951, the work of Bianchi was introduced in cosmology by Abraham Taub 
into his article "Empty Spacetimes Admitting has Three-Parameter Group 
of Motions" 112]. Spatially homogeneous spacetime have a time dependant 
space geometry which is thus a homogeneous 3-geometry. Thus, spacetime 
has a group of isometry with r dimensions acting on a family of hypersurface 
with r = 3 (simply transitive action), r = 4(locally rotationally symmetric) 
or r = 6 (isotropic models). The case r = 3 became known under the name 
of Bianchi cosmologies after the article of Taub. During one decade, the 
Bianchi models felt into forgets until the rebirth of General Relativity at 
the beginning of the sixties. O. Heckmann and E. Schiicking made them 
appear again in 1958 in their work "Gravitation, an Introduction to Current 
Research" Pni- Then it was with the turn of the Russian school of Lifshitz 
and Khalatnikov, joined later by Belinsky (BKL), through their study on the 
chaotic approach of the initial singularity which inspired Misner in the USA 
and later Hawking and Ellis in the United Kingdom. The Bianchi classifi- 
cation was reexamined by C.G. Behr in a work not published but reported 
in in 1968. Finally an essential contribution on the Bianchi models was 
made by Ryan, a student of Misner, and summarized through its book "Ho- 
mogenous Relativistic Cosmologies " pl| . 

Let us explain the physical reasons which, in our opinion, justify the study 
of these models. The Universe such as we observe it today is very well de- 
scribed by the isotropic and homogeneous cosmological models of Friedman- 
Lemaitre- Robertson- Walker (FLRW). It was shown by the observations of 
the cosmic microwave background (CMB) by the satellites COBE and WMAP[T9]. 
However nothing enables us to extrapolate these properties of isotropy and 
homogeneity at early times before the radiation/matter decoupling. The 
question thus arises of knowing why the Universe has these properties whereas 
there is an infinity of cosmological models not having them. Primarily, two 
answers are distinguished: 

• It could exist a quantum principle which selects, among the set of the 
possible models, the FLRW ones as being the most probable. This 
answer rests on the development of a quantum theory of the initial 

*It also exists the axisymmetric anisotropic model of Kantowski and Sachs 
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conditions. 



• the early Universe would be inhomogeneous and anisotropic but its 
evolution would lead it to a homogeneous and isotropic state (asymp- 
totic or temporary) corresponding to FLRW models. 

It is this second point of view which we will adopt by considering that the 
Universe is initially anisotropic and asymptotically becomes isotropic. We 
will keep the assumption of homogeneity because the inhomogeneous cos- 
mological models are not classified due to their lack of symmetry. Moreover 
we will consider that this state is reached asymptotically and not transitory. 
Indeed, the observations show that our Universe must be isotropic since at 
least its first million years and it is thus reasonable to suppose that once 
this state reached, it is stable. To consider that the Universe is initially 
anisotropic thus makes it possible to explain the isotropisation driving pro- 
cess rather than to consider this state in an ad hoc way. Another advantage 
of the anisotropic models resides in their approach of the singularity in Gen- 
eral relativity which is oscillatory and chaotic for the most general of them. 
It would be shared, according to the BKL conjecture, by the inhomogeneous 
and anisotropic modelsjSOl on the contrary to the FRLW models whose sin- 
gularity approach is monotonous. 

In this chapter, we will be interested by the isotropisation process for the 
Bianchi class A models with the aim to constrain the scalar-tensor theories 
so that they led the Universe to a stable isotropic state. For that we will use 
a method inspired from Wainwright, Ellis et al book "Dynamical Systems 
in Cosmology" |2i]. In this book, the dynamical evolution of homogeneous 
cosmologies without scalar field is described using first order equations sys- 
tems obtained thanks to the orthonormal formalism. These equations are 
then studied from the point of view of their equilibrium points. We will 
proceed in the same way but we will use the Arnowitt-Deser-Misner(ADM) 
Hamiltonian formalism and will seek the equilibrium points representing a 
stable isotropic state for the Universe. 

Most of the results which will follow were published in [22j 1221 Ej- This 
chapter is thus for us an occasion to gather them and to improve them con- 
siderably. Indeed, since these first papers on isotropisation, some significant 
progress which we will report here, were made. They consist of: 

• a classification of the Universe stable isotropic states in three classes. 

• a clarification of our results range thanks to a discussion on their sta- 
bility. 

• the use of numerical methods allowing to validate our calculations and 
to represent them graphically. 



5 



• some new physical interpretations of our calculations and results. 

The plan of this work is as follows. In the section 121 we mathematically 
describe the geometrical and physical frameworks of this work by writing 
the metric of the Bianchi models, the action and ADM Hamiltonian of the 
scalar-tensor theories. In the sections El and 01 we study the isotropisation 
of the Bianchi class A models. In the sectional we analyse the quintessent 
asymptotical nature of the scalar field. We conclude in the section 

2 Background 

This section mathematically describes the geometrical and physical frame- 
work of this work, namely the Bianchi models metric, the action and ADM 
Hamiltonian of the scalar-tensor theories and the method which we will sys- 
tematically use to study the isotropisation process. 

2.1 Metric 

The metric that we will use reflects the 3-1-1 decomposition of spacetime, 
necessary to obtain the ADM Hamiltonian. It is written: 

ds'^ = -(7V2 - NiN')dn'^ + 2Nidnuj' + RlgijJJ (1) 

being the lapse function, Ni the shift functions and the 1-forms gen- 
erating the Bianchi homogeneous spaces. For example for the spatially flat 
Bianchi type / model, we have uj^ = dx, a;^ = dy and = dz. We will 
choose a diagonal metric such as Aj = and the relation between the proper 
time t and the variable will be then: 

dt = -Ndn (2) 

2.2 Action 

The general form of the action of a minimally coupled and massive scalar- 
tensor theory with a perfect fluid is written: 

S = (Wtt)-^ j[R- (3/2 + c^)0'^0,^(A-2 - U\^d\ + S^iQij) (3) 

where the following functions are defined: 

• i? is the curvature scalar. 

• (/> is the scalar field. Here, it is minimally coupled because there is no 
coupling between (f) and the curvature and thus the gravitation constant 
does not vary with respect to this field. 
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• is the Brans-Dicke function describing the couphng of the scalar 
field with the metric. 

• U is the potential describing the coupling of the scalar field with itself. 
When the potential is not zero, the scalar field is known as massive 
otherwise it is massless. 

• Smigij) is the action representing the presence of a non tilted perfect 
fluid with an equation of state p = (7 — l)p, where p is its pressure, p 
its density and 7 its barotropic index such as 7 G [1, 2]. The case 7 = 1 
represents a dust fluid and 7 = 4/3, a radiative fluid. Note that Sm 
does not depend on (p. If it were the case, the matter will not follow 
the spacetime geodesies and it is likely that such a deviation would 
have been observed. 

The conservation of the perfect fluid energy-momentum tensor allows to get: 

where V = e~^^ is the Universe 3-volume. One will also define the pressure 
and the density of the scalar field like: 

^^ = 2^^^ +2^ (4) 

P^ = -2-^^ "2^ 
where a prime indicates a derivative with respect to the proper time t. 

2.3 Hamiltonian 

In this subsection we look for the ADM Hamiltonian H corresponding to 
the action lO). We follow the procedure described in PSj and [2n|- First, we 
rewrite the action as: 

S = (167r)"i j (ff^'^ + n<^^ - NC"^ - Ni&)d*x (6) 

The liij and II,^ are respectively the conjugate momenta of the metric func- 
tions and scalar field. The lapse and shift functions here play the role of 
Lagrange multipliers. By varying ^ with respect to and Aj, we get the 
constraints = and = with: 

Hj) 



i3)gU{cf>) + 



1 fe3(7-2)f7 
247r2 
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(8) 



where the "^^)" stands for the quantities defined on an hypersurface and 5 is 
a positive constant. We then rewrite the metric functions Qij as 



It means that stands for the isotropic part of the metric and I3ij for the 
anisotropic parts. Then, using the Misner parameterization |27j : 



p^fc = 27rnl - 2/3^5^ni 

Qpij = diag{p+ + V3p^,p+ - \/3p„, -2p+) 
I3ij = diag{l3+ + ^/3_ , [3+ - V3/3_ , -213+) 
we rewrite the action under the form 



(9) 
(10) 
(11) 



S = j p+df5+ +p-df3- + 



Hdn 



with p0 = vrll^ and H = 27rn^, the ADM Hamiltonian. Using the expression 
Q and the constraint C*^ = 0, we finally get the form of H: 



2 w!,2 



H'^=pI+pI + U^P^ + 247r2ii6e-6f^ C/ + 1/(17, /3+ , /3_ ) + 6e^^^'^^^ (12) 

where V{0,, (3+, I3-) stands for the curvature terms specifying each Bianchi 
model and given in table ^ 



Bianchi type 


m/3+,/3-) 


/ 





// 


127r2i?4e4(-Q+/3++V3/3-) 


VIo, VIIo 


247r2i?4e-4n+4/3+(pQgj^4^^_ ^ 


VIII, IX 


2A7r^R^e-^^[e^f^+icosh4:V3l3- - 1) + 
l/2e-8^+ ± 2e-2/^+ cosh2V3/3_] 



Table 1: Form of the curvature potential V{Q, P+, (3-) for each Bianchi model. 



2.4 Method 

The method that we will use to study the Bianchi models isotropisation 
process will be as follows: 

• Prom ADM Hamiltonian (|12p . we determine the Hamiltonian first order 
field equations. 
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• In order to use the dynamical systems methodspT], we rewrite these 
equations using bounded variables. 

• We then look for the equilibrium points representing the stable isotropic 
states. We deduce some constraints on the functions uj and U defining 
the properties of the scalar field and we calculate the exact solutions 
corresponding to the equilibrium points. 

• We apply our results to some scalar-tensor theories usually studied in 
the literature, which will enable us to test their validity. 



3 The flat Bianchi type / model 

The Bianchi type / model is the only one to be spatially flat. It thus contains 
the flat FLRW solutions. In the first subsection, we will study the isotropi- 
sation process without a perfect fluid and then, in a second subsection, we 
will consider a perfect fluid. We proceeded so that these two subsections can 
be approached separately: in each of them we wrote the field equations and 
discussed our results. 



3.1 Vacuum case 
3.1.1 Field equations 

The ADM Hamiltonian of the Bianchi type / model when no perfect fluid is 
present but with a minimally coupled and massive scalar field is written: 

2 i 2 

= + p2_ ^ u-M^ + 247r2i?6e-6f^^ (13) 
It comes then for the hamiltonian equations: 

■ dH 12(^2 



(3 + 2uj)H 



(15) 



- - 9^ - -^^3 + 2u;)F + ^^3 + 2^)2i7 " ^'^T- ^^^^ 
^-Mi'dn- "^^^ ^'^T ^^^^ 
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a dot meaning a derivative with respect to fJ. Having chosen Ni = and 
using the fact that d^/dQ = -l/2n^iV(2Sl, one deduces the form of the 
lapse function: 

iV=l?4^ (19) 
H 

The lapse function will have the same form in the presence of a perfect fluid 
or of curvature. In order to find the isotropic equilibrium points of this 
equations system, it is necessary to partially rewrite it using some bounded 
variables. The form of the ADM Hamiltonian suggests us defining: 

X = H-^ (20) 

y = 7rRlVe-^^UH~^ (21) 
z = p^(j){3 + 2ujy^^'^H-^ (22) 
These new variables have all a physical interpretation: 

• the variable is proportional to the shear parameter S introduced 
into (21]. 

• the variable is proportional to {p^p — p^) / [dO./ dt)"^ , (d^l/dt)'^ repre- 
senting the Hubble function when the Universe is isotropic. 

• the variable is proportional to {p^ + P(p)/{dQ/dt)'^ . 

• Defining the scalar field density parameter 0,^ like the ratio of its energy 
density with the square of the Hubble function, we deduce from the 
last two points that 

— is a linear combination of and z'^ . 

— is proportional to and when the scalar field is quintessent. 

One can check that these variables are bounded by rewriting the Hamiltonian 
constraint in the form of a sum of their squares: 

+ 24y^ + 122^ = 1 (23) 

with = p\ + p"^. So that they be real, we will consider that 3 + 2u; and 
U are positive functions. The first condition is necessary to the respect of 
the weak energy condition which states that +P0 > 0. The Hamiltonian 
equations can be rewritten according to these variables in the following first 
order differential equations system: 

X = 72y^x (24) 

y = y{6ez + 72y2 _ 3) (25) 
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i = 24!/2(3z _ If) (26) 

where £ is a function of the scalar field defined by: 

^ = (/,[/^C/-i(3 + 2a;)-i/2 

We thus reduced the seven Hamiltonian equations in the form of a system 
of four equations with four variables x, y, z and (/>, this last one not being 
necessarily bounded. This reduction of the equations number comes owing 
to the fact that the equations of Hamilton show that ■p± are some constants 
and that /?+ oc One thus eliminates three from the seven equations of 
Hamilton. Our next objective is then to find the equilibrium points corre- 
sponding to a stable isotropic state for the Universe and whose properties 
are defined in the following section 



3.1.2 Definition of a stable isotropic state 

Isotropy is defined in the article of Collins and Hawking[2Hj when the proper 
time diverges using the four following properties: 

• — > — oo 

This condition says to us that the Universe is forever expanding when 
isotropisation occurs. Considering that no period of contraction was 
observed since radiation-matter decoupling and that our Universe is 
strongly isotropic, this appears justified. 

• Let be the energy-momemtum tensor: T^^ > and — > 

^5fj represents a mean velocity of the matter compared to surfaces of 
homogeneity. If this quantity did not tend to zero, the Universe would 
not appear homogeneous and isotropic. 

• Let be aij = {de^ /dt)k(i{e~'^)j)k and cj^ = o-jjcr*-?': 0. 

This condition says that the anisotropy measured locally through the 
constant of Hubble Hq tends to zero. Indeed, when we measure the 
Hubble constant, we evaluate the quantity Hi = ^j^/gu = djda/dt — 
dQ/dt. So that it appears to be the same one in all the directions, it 
is thus necessary that dPa/dt « dVi/dt. 

• P tends to a constant Po 

This condition is justified by the fact that the anisotropy measured 
in the CMB is to some extent a measurement of the change of the 
matrix f3 between time when radiation was emitted and time when it 
was observed. If (3 did not tend to a constant, one would expect large 
quantities of anisotropics in some directions. 
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Within the framework of the Bianchi type / model, we have that dj3±/dt = 
—N~^(3± oc e^^. Hence, so that j3± tend asymptotically to a constant in ac- 
cordance with the fourth point, it is necessary that Vt, — > — oo in agreement 
with the first one and then, d/3±/dt vanish. However, generally, when the 
derivative of a function tends to zero in t — > +00, this necessarily does not 
imply that the function tends to a constant. It is for example the case of 
the logarithm Int. Hence the last point in Collins and Hawking definition 
implies that the isotropy appears relatively quickly. Note that the vanishing 
of dP±/dt is necessary so that the Universe be isotropic otherwise. Hi — H2 
as instance will be different from zero which is not observed. This is an 
important point because in presence of curvature, we will see that the limit 
f2 — >■ — 00 is not sufficient so that dp±/dt vanish whereas here it is essential 
and is a better justification than the convergence of f3± to a constant. Con- 
cerning the third point, it means that the shear parameter proportional to 
X (X $± = tends to zero. Consequently, the stable isotropic equilib- 

rium points that we are looking for will be such as: 

O — > —00 

x^O 

By examining the equations system H24II26|1 . one notes that there are three 
types of isotropic equilibrium states corresponding to these two limits and 
making it possible to define three isotropisation classes in the following way: 

• Class 1: all the variables {x,y,z) reach an isotropic equilibrium state 
with y 7^ 0. It is the class which seems to correspond to the most 
studied scalar-tensor theories in the literature. 

• Class 2: all the variables {x,y,z) reach an isotropic equilibrium state 
with y = 0. In this case, it is generally not possible to determine the 
asymptotic behavior of x with the approach of the isotropy because 
it depends on the unknown way in which y tends to zero whereas it 
is X which determine the common asymptotic behavior of the metric 
functions at the isotropisation time. The class 2 isotropisation was 
observed numerically for a non minimally scalar field in |29j . 

• Class 3: x tends to equilibrium but not necessarily the other variables. 
As they must be bounded when Q — > —00, that means that they have 
to oscillate in such a way that their first derivative with respect to Q 
oscillate around zero. The class 3 isotropisation was observed numeri- 
cally for complex scalar fields in (HH] , 

In this work our attention will be concerned with the study of class 1 isotropi- 
sation. 
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3.1.3 Study of the isotropic equilibrium states 

This section breaks up into three parts. In the first one, we calculate the 
location of the isotropic equilibrium points and in the second one we describe 
the assumptions which will enable us to derive the asymptotic behaviors of 
the metric functions and potential in the vicinity of the isotropy in the third 
section. 

equilibrium points 

The equilibrium points of the system (|24II26|) such as a; = and y ^ 
are given by: 

(x, y, z) = (0, ±(3 - ^2)1/^(6^2), ^/6) (28) 

In order to reach equilibrium, it is thus necessary that tends to a constant 
smaller than three. Linearising the equation (j24p in the vicinity of equilib- 
rium, it is found that asymptotically when — > — oo, the variable x behaves 
like: 

^^^oe3C-/£2dn (29) 
Assumptions and stability of the asymptotic behaviors 

Before further going, this first calculation (|29|l of an asymptotic behavior 
offers the opportunity to speak about the stability of our results. They can 
be separate in three types: 

1. The location of the isotropic equilibrium points. 

2. The conditions necessary to their existence. 

3. The asymptotic solutions associated to the equilibrium points. 

However the asymptotical solutions and some of the necessary conditions will 
be determined by calculating the exact solutions for each equilibrium point 
and they will be correct only if on one hand i and in the other hand the vari- 
ables (y, z) tend sufficiently quickly to their equilibrium values. Otherwise, 
they will be different and we are going to explain why. 

With regard to £, in order to obtain some of our results in closed form and 
comparable between them, we will generally make the following assumption 
that we will call "variability assumption of £": 

• When near equilibrium £^ tends to a constant (.q, vanishing or not, 
with a variation 5f such as ^ ll + 6f, J fdQ + const (If 

^0 = 0, the integral thus tends to a constant). 

We will see in the third part of this section that it is possible to asymp- 
totically determine the behavior of 4>{^) and thus of Consequently, it 
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is possible to avoid any approximation with respect to I as calculation l(29|l 
shows it above: whatever the speed to which I tends to its equilibrium value, 
the presence of the / I'^dQ term makes it possible to take into account the 
variation (5£^ of £^ in its neighborhood. In order to show that one can math- 
ematically free oneself from this assumption, all the results of the sectional 
will be expressed by taking into account the integral of l"^ and then the vari- 
ability assumption of I. Moreover, we will apply these results to the cases of 
two scalar-tensor theories, respectively in agreement and disagreement with 
this assumption. 

With regard to y and z, we will make the same type of assumptions as 
for I by considering that their variations 6y and 6z when we approach the 
equilibrium are sufficiently small to be negligible. For example, when we 
calculate the asymptotical behavior (|29|l for x, we take into account the way 
in which £ approach its asymptotic value since we do not make the variability 
assumption of I but not that of y which we simply replaced by its equilibrium 
value in the field equations without considering 5y. This problem can not 
be solved as "easily" as the one in connection with (.. It is possible that a 
perturbative study can bring some answers but it is not guaranteed because 
it could strongly depend on the specification of the forms of uj and U as 
functions of the scalar field. 

To summarize, all our results implying the calculation of an asymptotic 
approach of a quantity in the vicinity of equilibrium will be valid when the 
Universe reaches sufficiently quickly the isotropic state. The assumption 
on i can be raised while not making the variability assumption but that 
seems more difficult for the variables y and z. Hence, we will systematically 
assume that these variables sufficiently quickly approach their equilibrium 
values. These assumptions will be also valid for the variables k and w that 
we will define later, respectively associated with the presence of a perfect 
fiuid and curvature. 

Asymptotical behaviors 

Applying the variability assumption of (. to ^^I^ we have that when (. tends to 
a non vanishing constant, x —>■ e^^~^ and to e^^ otherwise. This variable 
thus disappears when fl — > — oo and the reality condition of the equilibrium 
points is respected. 

At the same time, the equation for x shows that x is a monotonous 
function of 0: when x is initially positive (negative), it is asymptotically 
increasing (decreasing), x is thus also of constant sign. Using the expression 
ifT!!]! for the lapse function N and owing to the fact that dt = —NdQ, we 
deduce that Q{t) is a decreasing (increasing) function of the proper time 
t when a; , or in an equivalent way the Hamiltonian, is initially positive 
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(negative). Consequently, an initially positive Hamiltonian is a necessary 
initial condition so that isotropy when Q — > — oo occurs at late times. Fi- 
nally, a last remark about the monotonous functions. We can calculate that 
dgij/dQ = — 2e~^^"'"^*-' (1 — Taking into account what we said on the 
monotony of x and the expression for the derivatives of Pij with respect to 
O, it comes that the Pij are monotonous functions of proper time t and that 
consequently each metric function can not have more than one extremum 
during time evolution. We had shown this point in a different way in by 
using the ADM Hamiltonian formalism. 

To determine (/'(fi), we use the equation l(27|l for cp which is written asymp- 
totically: 

^ C/(3 + 2lo) 

In this last expression the variability assumption of i is not made but on the 
other hand we neglect the variation 6z of z in the isotropy neighborhood. 
It is the asymptotic form of the solution of this equation which will give us 
the asymptotic behavior of (p according to O. One will thus deduce £{Q) 
and [/(O) which are two unknown functions of the scalar field. In particular, 
to know £(0) will make it possible to check the necessary conditions for 
isotropy when O — > — oo, the variability assumption of £ and to calculate 
the asymptotic form of the metric functions Q{t) and the potential U{t). 
Indeed, using on the one hand the asymptotic behavior of x and the relation 
dt = —NdQ and on the other hand the definition of y , it is found that with 
the approach of the isotropy Q and U behave respectively like: 

dt = -UnRlxoe-f^^'^^dQ. 

and 

and then, if we take into account of the variability assumption of i 

dt = -UnRlxoe-^^^dn 

and 

U = e^^'^ 

This assumption enables us to calculate that when i'^ tends to a non vanishing 
constant, the metric functions tend to t^ and the potential to On the 
other hand, when £^ tends to vanish, the Universe tends to a De Sitter 
model and the potential to a cosmological constant. If the assumption is not 
checked, it is always possible to determine these asymptotic behaviors but 
with some quadratures. 
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3.1.4 Discussion and applications 

The results of section IXTl are concerned by a minimally coupled and massive 
scalar-tensor theory without a perfect fluid. It is the simplest of the theories 
that we will consider and the method used above will serve us as a guide for 
the next sections. We have assumed that U and 3 + 2a; are positive functions 
and a class 1 isotropisation. When it is supposed that the function £ and 
the variables y and z tend sufficiently quickly to their equilibrium values, we 
have the following results: 

Let us consider a minimally coupled and massive scalar-tensor theory and 
the function I defined by i = jjj^^^/yTJ- The asymptotic behavior of the 
scalar field to the approach of the isotropy is given by the form of the solu- 
tion when 0, — > — oo of the differential equation <j) = '^ u(3+2uj) ■ ^ necessary 
condition for class 1 isotropisation is that £^ < 3. If £ tends to a non vanish- 
ing constant, the metric functions tend to and the potential disappears 
like If £ tends to vanish, the Universe tends to a De Sitter model and 
the potential to a constant. 

The variability assumption of I can be raised and the results are expressed 
then using the integral of £^ as shown above. They are in agreement with the 
"Cosmic No Hair theorem" of Wald^^ which says that homogeneous and 
initially expanding models with a positive cosmological constant (except the 
Bianchi type IX model) and an energy-momemtum tensor satisfying the 
strong and dominant energy conditions, tend to an isotropic De-Siter model 
for which the expansion is exponential. Here, when a cosmological constant is 
considered or when — > such as the variability assumption of I is checked, 
the Universe, when it isotropises, tends to a De Sitter model and the po- 
tential to a constant. On the other hand, when I tends to vanish but the 
variability assumption of i is not checked, the potential does not tend any 
more to a constant and the Universe does not approach a De Sitter model. 

As an application, we will examine the cases of the scalar-tensor theories 
defined by the Brans-Dicke coupling function 

and the forms of potentials 
and 

u = (r 

It is pointed out that some of our results represent necessary conditions 
and that consequently, when in the applications below we will speak about 
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isotropisation, it is always by keeping in mind that they could not be suffi- 
cient. 

The exponential potential of (f) was largely used in the literature. Isotropi- 
sation of Bianchi models with this potential was already studied in jSHl and 
thus will enable us to test our results. It was shown that all the Bianchi 
models (except the Bianchi type IX model when it contracts) isotropised 
at late times when rm? < 2. If m = , the Universe tends to a De Sitter 
model because the potential is a constant and, if not, it is expanding such as 
t^™ . linn? > 2 , the Bianchi type /, V , VII and IX models may 
isotropise at late times. Now, using our results, we get that asymptotically: 

(j) — > mVL 

The necessary condition for class 1 isotropisation is then m? < 6 and the 
asymptotic behaviors of the metric functions are well in agreement with 
what was predicted in [3S1- The difference between the results of this last 
paper and ours is the nature of the range of m authorizing the isotropisation 
since we find an upper limit. The figure ^ illustrates the convergence of the 
variables x, y and z to their equilibrium values for m = — 1. When rm? > 6 



y/(iiE0'3) 




Figure 1: Evolution of the variables x, y and z when (^+^^)^''^ = U = e'"'*, i?g = 

(V247r)-l and m = -1 with the in itial values ix,y,z,4,) = (0.87,0.25,-0.12,0.14). £ = ~1/V2, 
X tends to 0, ?/(7rRg)-l to VS - P/{6ttRIV2) = 0.91 and z to i/G = 0.12 in agreement with the 
equilibrium points expression. 

, the class 1 isotropisation is not possible any more because y equilibrium 
value would be complex. A numeric simulation of this case is shown on the 
figure El when m = —3.2, the Universe always tends to an equilibrium state 
which is now anisotropic because x tends to a non vanishing constant and 
thus the functions P± describing the isotropy diverge. 

Now let us examine the case of a power law potential of the scalar field. 
Then we get for i: 




and, if a class 1 isotropisation occurs, the scalar field behaves as 

We must thus have m < so that the scalar field is real and we deduce 
that £^ tends to vanish like m(40)~^. Consequently, / I'^dQ does not tend 
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Figure 2: Evolution of the variables x, y and z when (^+2^)'^^ = ^ = ^ ^3 ^ 

(y247r)-i and m = -3.2 with the initial values {x,y,z,,j>) = (0.87,0.25,-0.12,0.14). Universe 
does not isotropise: the system tends to an anisotropic equilibrium point such as the functions 
f3± describing the anisotropy diverge. 



to a constant when Q — > — oo but diverges like ^ ln(— fi) and we must 
take account of this integral in our results: the variability assumption of 
£ is not checked here. Raising this assumption, one finds then that the 
potential tends to vanish like (— fi)™/^ and the metric functions tend to 



exp 



( ^ S — t)*-'" . So that this quantity diverges positively, it is thus 



necessary that m < 4 what is always checked since m < 0. This case is 
illustrated on the figure 01 where it is seen very clearly that convergence 
of the variables y and z to their equilibrium values is slow (compare with 
the figure ^ of the previous application) . This means that the Universe 
approaches "slowly" its isotropic state and one could then think that, in 
addition to raising the variability assumption of £, the variations Sy and 
6z of the variables y and z near equilibrium about which we speak in the 
previous sub-section should also be taken into account. However, it seems 
that these last corrections are not necessary. This can for instance be checked 
by comparing the asymptotic evolution of z corresponding theoretically to 
z ^ i/6 with the numerical integration of the figure 01 for the large Q values. 
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Figure 3: Evolution of the variables x, y and z when (3+^^)'^^ = ^ [/ = ^3 ^ [^t,)-^ 
and m = -1 with the initial values {x,y,z,(f>) = (0.87,0.25,-0.12,0.14). x tends to 0, 1/(71 i?g)-i 
to 1 and 2 to in agreement with the equilibrium points expression. Let us notice that the 
variables y and z tend much less quickly to their equilibrium values than on the figure |2] This is 
due to the "slowness" of the convergence of £ to zero which then reflects on the variation of these 
variables. 

When m > , a class 1 isotropisation does not seem possible any more 
because then the scalar field would be complex. Numerical integrations show 
that the initially positive scalar field decreases to zero for a finite O value. 
Then they fail for negative value, maybe indicating a singularity occurrence. 
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3.2 Perfect fluid case 



The method is the same one as without a perfect fluid but an additional 
term appears in the field equations|23] due to its presence. Its equation of 
state is p = (7 — l)p with 7 G [1,2]. We will also consider that the pres- 
sure of the non tilted perfect fluid is isotropic. It is a simplifying assump- 
tion whose consequence for the Bianchi type / model could be to increase 
the anisotropics decay(in V~^), thus preventing any observational detec- 
tion, whereas the presence of an anisotropic pressure slow down it. In this 
last case, the anisotropy could then be detected via the ratio 5T/T of the 
CMB which depends on the shear parameter at the time of last scattering 
surface [H]. 



3.2.1 Field equations 

This time the ADM Hamiltonian is written: 

H'=pI+ pi + 12^^ + 2^7:^ Rle-^^'U + de^^^-^^"" (30) 
+ 3 + 2a; " ^ ' 

where (5 is a positive constant. Compared to the previous section, one thus 
sees appearing the term 5e^'^'^~'^^^ due to the presence of the perfect fluid. 
The Hamiltonian equations become: 

■ dH 12(/)2p, 



(3 + 2uj)H 



dH 



(?>^2<jS)H (3 + 2w)2i/ 'J H 

^JT f)TT -&Q.TJ 3(7-2)n 

In order to rewrite these equations, we will use the same bounded variables 
X, y and z that previously defined to which we will add a fourth variable: 

^2 ^ ^g3{7-2)n^-2 



related to the presence of the perfect fluid. This variable is in fact pro- 
portional to the density parameter of the perfect fluid, one of the main 
cosmological parameters usually noted f^m- This can be shown by checking 
that k"^ oc y~'^' /{^Y where ^ is the Hubble function when the Universe 
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isotropises. k is not independent of the other variables and can be rewritten 
hke: 

= fe2e3(^-2)f^ (32) 
= 6y^U~W-^ (33) 
Then, the Hamiltonian constraint becomes: 

+ R^y^ + 12^2 + = 1 (34) 

and the Hamiltonian equations are: 

X = 72y^x - 3/2(7 - 2)k'^x (35) 

y = yiGiz + 72y2 _ 3) - 3/2(7 - 2)k^y (36) 
f 

z = 2V(3^ - -) - 3/2(7 - 2)k'^z (37) 

(/) = 12z- (38) 

^ (3 + 2a;)V2 ^ ' 

with always £ = (pU^U'^i^ + 200)-^/"^. 
3.2.2 Study of the isotropic states 

One distinguishes two types of equilibrium states according to whether k , 
i.e. the density parameter of the perfect fluid, tends to vanish or to a non 
vanishing constant. 

k tends to vanish 

As y does not tend to vanish because we consider a class 1 isotropisation, 
we deduce from the form l(33)l of k that U » V^^. The equilibrium points 
are the same ones as in the absence of perfect fluid and thus we find again 
the reality condition < 3. The asymptotic behavior of x near equilibrium 
is obtained starting from the equation Ip^ : 

From this expression and the definition of y, we then deduce for the potential: 



While making use of the asymptotic behavior of x and the definition (j32|l of 
k, it comes: 
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By deriving this expression, one obtains the differential equation 
2kk = _ _ 2)/c2 ^ 3J ^2 



whose exact solution is: 

^2 ^ ^ 

/cq being an integration constant. All these results were obtained without 
applying the variability assumption of I. So now we take it into account and 
we obtain: 

It thus follows that /c — > when Q. — —00 if < ^ < 3. Consequently the 
isotropisation will occur when A; ^ only if £^ < ^ which is a smaller range 
than when no perfect fluid is present(then £^ < 3). However, the asymptotic 
behaviors of the metric functions and the potential remain unchanged but if 
k sufHciently quickly does not approach its equilibrium value. 
If the variability assumption is not valid, again one must take into account 
the integral of The range of £^ allowing the isotropisation will be always 
at least such as £^ < 3 because this inequality is independent of any ap- 
proximation but it will be modified (or not) differently by the fact that we 
consider the limit /c ^ 0. Moreover the asymptotic behaviors of the metric 
functions and potential will be different from what they are when there is 
no perfect fluid in spite of this disappearance of k. 



k tends to a non vanishing constant 



The isotropic equilibrium points are not more the same ones and thus ei- 
ther the asymptotic behaviors of the metric functions and the potential. For 
the first ones, we find: 



after having deduced from the constraint that 

A;2 = 1 - 37(2^2)-! 

The equilibrium points will be real if 7 is a positive constant smaller than 
2, in agreement with the range of variation of 7 which we specified, that 
is 7 G [1,2]. The variable k will be real and the other variables will reach 
equilibrium for a not vanishing value of y, thus respecting the definition of 
class 1 isotropy, if i'^ tends to a constant larger than This necessary 
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condition for isotropy is independent of any approximation. Linearising the 
differential equation for x , we find that asymptotically: 

and that the metric functions tend to 

From the definition of y, we deduce that the potential tends to vanish as 
what is confirmed by the form of k which shows that asymptotically 

U oc 

in agreement with the asymptotic expressions of the potential and the metric 
functions with respect to the proper time t. This last expression makes it 
possible to determine the scalar field asymptotic form according to the 
one of the potential. Let us note that all these asymptotic behaviors are 
independent of the variability assumption of 

3.2.3 Discussion and applications 

Let us summarize our results obtained in the presence of a perfect fluid. For 
that, we will state them according to the density parameter of the perfect 
fluid which is proportional to k. When it tends to vanish, we will make the 
variability assumption of I otherwise it is useless as shown above. It comes: 

Isotropisation with ^Im 0: 

The results are the same ones as without a perfect fluid but the range of £^ 
allowing the isotropisation is reduced to i'^ < Moreover, at the time of 
the isotropisation the scalar field potential is asymptotically larger than the 
perfect fluid energy density. 

When the variability assumption of i is not made, the things are not so 
simple and the results depend completely on the way in which i"^ approaches 
equilibrium. One can however always calculate them by using the expres- 
sions given in the preceding sections and depending on / ^^dO. 

When k tends to a non vanishing constant, the equilibrium state is different 
and we find that: 

Isotropisation with Qm -h O-' 

Let us consider a minimally coupled and massive scalar-tensor theory and the 
quantity i defined by i = jj^^^y^TTJ- The scalar field asymptotic behavior at 
the approach of isotropy can be deduced owing to the fact that U oc : the 
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potential of the scalar field is proportional to the perfect fluid energy density. 
A necessary condition for a class 1 isotropisation will he that i"^ > i'^ 

is bounded and < 7 < 2. Then the potential disappears like and the 

2 

metric functions tend to t^t . 

This last result shows that the scalar-tensor theory tends to General Rel- 
ativity with a perfect fluid in agreement with [SS]- Then, the presence of the 
scalar field have no effect on the asymptotic evolution of the metric func- 
tions. Thus for a fluid of dust such as 7 = 1, the Universe tends to that of 
Einstein-De Sitter with — > t^/^ and for a radiative fluid such as 7 = 4/3, 
to a Tolman Universe with t^/"^ . Also let us notice that the ranges of 

l'^ allowing isotropisation when A; ^ and k -/^ are complementary. 

Now we examine our results in the light of the two scalar-tensor theories 
which we had considered in the absence of a perfect fluid in the section 

ura 

Let us start by considering an exponential potential. The only difference 
with the results obtained without a perfect fluid is the necessary condition 
allowing the isotropisation. When m? < 37 , isotropisation should occur 
with Qrn — *■ 0, and the metric functions tend to t'^"^ . This case is illus- 
trated on the figure @] for m = 1. When > 37 , isotropisation should 

_2_ 

occur with A 0) ^iid the metric functions tend to t^T . This case is illus- 
trated on the figure ini for m = 4. These results are in agreements with those 
found in jHEl for the FLRW models. However, in this last paper, a stable 
solution of trackers type such as Qm — *■ had also been found when > 6. 
Here, we do not find it since it does not allow for isotropy. 
With regard to the power law potential when 0, we know that the 



Figure 4: Evolution of the variables x , y, z and k in the presence of a fluid of dust (7 = 1) 
when (3+2^)^/^ = = e™'/', r3 = (v^tt)'! and m = 1 with the initial values {x, y, z, (f>) = 

(0.62,0.25, —0.12, —0.14). All occurs as in the vacuum case with k (or equivalently Qm) tending 
to zero. 

variability assumption of £ is not verified since £ tends to zero but that the 
integral of its square diverges like ^ln(— fi) when Q — > — oo. Moreover, we 
showed that m must be negative so that the scalar field be real. Taking into 
account these elements, the calculation of k'^ and its integral gives us then: 
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Figure 5: Evolution of the variables x, y, z and k in the presence of a fluid of dust (7 = 1) 
when (3+2^)1/^ = 72, (7 = e""l' , = (v^tt)"! and m = 4 with the initial values (x, y, z, if)) = 
(0.62,0.25, —0.12, —0.14). y, z and k oscillate to their equilibrium values which are respectively 
-0.0625, -0.088 and 0.90 when £ = -4/^2. 



and 

J k'^dn oc r(l - m/2, -37O) ^ 

when Q — > — 00, r being the Euler function. Consequently, these two quan- 
tities also tend to zero when i7 — > — 00 without one needing additional con- 
dition. Simulations with m = — 1 and a dust fluid give identical results to 
the figures El 

When Qrn tends to a non vanishing constant, the scalar field behaves like 
(j) —> e rn and thus i disappears or diverge, preventing a class 1 isotropisa- 
tion. 



4 Bianchi class A models with curvature 

In this section, we will consider the presence of curvature by studying the 
Bianchi class A models of types //, VIq, VIIq, VIII and IX without a 
perfect fluid. This last model in particular, contains the solutions of the 
FLRW models with positive curvature. 



4.1 Field equations 



The ADM Hamiltonian for the Bianchi class A models with curvature is 
written: 

.2 J,2 



H'=pi+pi + 12^^^ + 2A7T'R^oe-'"'U + V{n,(3+,p^) (40) 

where /?+,/?_) is the curvature potential characterizing each Bianchi 
model and is given in tabled The Hamiltonian equations are then: 

dH_ 124>^p, 

dp^ - {?, + 2u)H ^^^^ 
dH dV 

= "a^ = -map: ^^'^ 
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- " 9^ - -^^3 + 2u;)F + ^^3 + 2^)2if " ^^^^ 

As for the Bianchi type / model, the isotropy means that dp±/dCl have to 
vanish. But from now on the conjugate moments of the P± functions are not 
any more some constants and thus df3±/dt = p±e^^: contrary to the Bianchi 
type / model, the limit — > — oo does not imply any more that d(3±/dt 
and thus that isotropy occurs for a forever expanding Universe. We have to 
show it. Let us suppose that isotropisation leads to a static Universe, i.e. 
such as O — > const, when the proper time diverges. Then, so that dj3±/dt 
disappear, it is necessary that p± 0. However the equations indicates 
that dp±/dt oc Consequently if O and f3± tend to some constants 

when t — > +00, for the Bianchi types II, VIq and VIII models, these deriva- 
tive tend to non vanishing constants and the conjugate momenta p± can not 
disappear and the isotropy to occur. On the other hand, the things are not 
so simple for Bianchi type VIIq and IX models because if (3± — > 0, it is the 
same for and one can say nothing on the asymptotic values of p± . We 
will further show for these models also, that the isotropy can emerge only 
for a diverging value of il. 

Consequently, for the Bianchi models with curvature, isotropisation can oc- 
cur only if: 

O — > iboo 

df5± 



dn 







We will see that the third limits are always satisfied when it is also the case 
for the two first ones. In what follows, the variability assumption of £^ will 
be systematically applied. We did not explore what occurs when that it is 
raised. Indeed, we will see that the results obtained for the isotropisation of 
models with curvature are similar to those obtained for a flat model. However 
the conditions to check to show that the isotropy is reached are much more 
numerous and the variability assumption can not be raised easily. Let us 
note however that this is technically feasible in the same way as for the 
Bianchi type / model. 

In order to describe the curvature of the Bianchi models, we will introduce 
new variables prefixed w and similar to the three Nt variables defined by 
arguments of symmetry of the structure constants in j2Tj and |37j . 

4.2 Bianchi type // model 

In order to rewrite the field equations, we use the following variables: 

x± = p±H-^ (46) 
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y = TTRlVUe-^^H-^ (47) 

z = p^0(3 + 2cu)-i/2^-i (48) 

w = 7ri?2e-2f^+2(/3++v^/5-)^-i (49) 

Only one w variable is sufficient to describe the curvature just as only one 
variable Ni was sufficient in j^. Then the condition — > 0, necessary 
to the isotropisation, writes x± — > 0. It will be the same one for all the 
Bianchi models for whom we will re-use the same variables x±, y and z. The 
Hamiltonian constraint and the field equations rewrite like: 

xl + xi+ 24y2 + 12^2 + 12^2 = 1 (50) 

x+ = 72y\+ + 2Aw'^x+ - (51) 

X- = 72y^x^ + 24w^x^ - IAVSw^ (52) 

y = y{6iz + 72y2 - 3 + 24^^) (53) 

z = y^(72z - 121) + 24w^z (54) 

w = 2w{x+ + VSx- + 12u;2 + 36^^ - 1) (55) 

The constraint shows that the variables (|46l49p are bounded. Moreover, we 
find the usual equation for the scalar field: 

■ 12(t> 



V3 + 2a; 



In order to determine the asymptotic behavior of the metric functions and 
potential, we will need to know the asymptotic behavior of the Hamiltonian. 
We thus rewrite the Hamiltonian equation for H in the form: 

H = -H{72y^ + 2W) (56) 

It shows that H is a monotonous function keeping its initial sign during its 
evolution. Consequently, one deduces from the lapse function that when H 
is initially positive (negative), O — > — oo corresponds to late (respectively 
early) times and vice versa when Q — > +00. 

With these equations, we can from now on calculate the equilibrium points 
of the equations (|51l55p . It comes: 



{y,w) = (0,0) 

{x+,x^,y,z,w) = (l,\/3,0,0,±2/2) 
{x+,x_,y,z,w) = {0,0,^^^,e/6,0) 
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(/, ^, — V|T^, V':>|T:j:g, =n 2(^;^+8) ' 2(^^+8)'^ 2(F+8) 

The first one is such as y = and can thus be ehminated because it does 
not correspond to a class 1 isotropisation. The second and the last ones do 
not allow to obtain x± — > 0. Consequently only the third equilibrium point 
corresponds to an isotropic state if i tends to a constant such as £^ < 3. 
It is similar to the one found for the flat Bianchi type / model. In order 
to find the asymptotical behavior of w, we linearize l|55p in the vicinity of 
equilibrium by neglecting the variables w and x± tending to zero. It comes: 

Linearising in the same way the equations H51I52|1 for x± , using the variability 
assumption of £^ and introducing this last expression for w , we obtain that 
x± behave as the sum of two terms e^*^^~^ and e^^~^ Since isotropy 
needs x± — > and £^ < 3, we deduce that it occurs only when £^ < 1 in 
O — > — oo. The special value £^ = 1 is not compatible with the isotropy. 
It may be deduced from our variability assumption of £^ which implies that 
if £^ — > 1, £^ — 1 disappears generally more quickly than But then, 

w would tend to a non vanishing constant which is incompatible with the 
expression of the equilibrium points. Finally we find that asymptotically 

The two limits £^ < 1 and Q —>■ — oo allow x± but also w to vanish: when 
the Universe isotropises, it is forever expanding. In order to know if our 
model may isotropise, it is also necessary to check that p±e^^ — > when 
O — > — oo. For that, we write p±/H like a function of x± and w and we 
use their asymptotic behaviors. One then calculates that p±/p± tend to 
the constant —(1 + ^^). Consequently, p±e^^ — > e*^^~^^)^ and disappear 
when VL diverges negatively and the necessary conditions for isotropy are 
respected. The asymptotic behaviors of the metric functions and potential 
are the same ones as for the Bianchi type / model and depend on the same 
manner of disappearance or not of the function £^ to the approach of isotropy. 
Concerning the 3-curvature, it tends to zero when ^ — oo, showing that 
the Universe becomes flat. 

4.3 Bianchi type VIq and VIIq models 
Now the variables which we will use are: 

x± = p±H~^ (57) 
y = TiRle-^^U^/'^H-^ (58) 
^ = p^</,(3 + 2J^)-V2//-i (59) 
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(60) 



The difference with the Bianchi type II model is that we need two variables 
Wi to describe the curvature, one of them {wj^) being the variable w previ- 
ously defined for this last model. This is in agreement with j2l] where two 
variables Ni are also necessary for these models. The Hamiltonian constraint 
is written: 



x_ 



+ xi + 2V + + 12(w;+ ± W-Y = 1 



and the field equations become 

x+ = 72y'^x+ + 24(x+ - l)(w_ ± w+)' 
X- = 72y^x^ + 24x_(it;_ ± w+ f + 24^/3{w'^_ 
y = y{6iz + 72y2 _ 3 + 2A{w- ± w+f 
i = y2(72z - m) + 2Az{w^ ± w+f 
+ + ^x^ + I2{w 



2w^ 
2w_ X 



± w+f + 36y2 



^/32;_ + 12('u;_ ± w+f + 36y2 



(61) 

(62) 
(63) 
(64) 
(65) 
(66) 

(67) 



Anew we express the Hamiltonian equation for H according to the variables 
y and w±: 



H = -H 



72y^ + 2A{w+±w^y 



(6J 



In the equations, the symbols it correspond respectively to the models of 
Bianchi type VIq and VIIq. For the first model, the constraint l(6T|l shows 
that the variables are bounded. It is not the case for the second one: be- 
cause of the minus sign, w+ and W- could diverge if the difference w+ — W- 
remains finished, thus respecting the constraint. We will show at the end of 
this section that this is in fact impossible. Supposing that all the variables 
(but 15^!) are bounded allow to derive that isotropisation is impossible for a 
bounded value of Q.. Indeed, if 17 — > const when the proper time t diverges, 
dVL/dt 0. But from the form of the lapse function and owing to the fact 
that dt = —NdQ , we deduce that H should also tend to zero. It comes 
then from the definition of the variables ■w± that they should diverge which 
is incompatible with the fact that they are bounded in the vicinity of the 
isotropy. Thus, the isotropisation can not lead the Universe to a static state 
and O diverges inevitably. 

We calculate the equilibrium points of the equations system (|62l67p . It 

comes: 

. (2/, «;+,«;_) = (0,0,0) 

• {x+,x^,y,w+,w-) = {1,0,0, w+,w+) 
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• {x+,x-,y,z,w+,w- 

• {x+,x-,y,z,w+,w- 

• {x+,x-,y,z,w+,w- 

• {x+,x-,y,z,w+,w- 

• (x+,0,0,0,±z/2) 

• {x+,x-,y,z,w+,w- 

• {x+,x-,y,z,w+,w- 



• {x+,x-,y,z,w+,w- 



• {x+,x-,y,z,w+,w- 

For the same reasons as for the Bianchi type // model, the only equilibrium 
points corresponding to an isotropic state are: 

(x+,x_,y, Z,w±) = {0,0,±VS- i^{6V2)-\i/6,0) 

It will be real if £^ tends to a constant smaller than 3. 

In the same way as for the Bianchi type II model, one can show that J7 is 
a monotonous function of the proper time whose divergence in — oo corre- 
sponds to the late times if the Hamiltonian is initially positive. We also got 
that the asymptotic behaviors of the functions x±, w± , p±e^^, and U 
are asymptotically the same, and thus that i"^ < 1 when we approach the 
isotropy. 

All these results was shown not by considering the individual behaviors of 
W-^- and W- but by considering that u>+ ± W- — > 0. As we deduce from 
this single limit that near the isotropic state w± — > 0, it follows that these 
variables are always bounded as stated at the beginning of this section, and 
in particular for the Bianchi type VIIq model. 

4.4 Bianchi type VIII and IX models 

We will use the following variables: 

x± = p±H~^ 

y = ^Jl3,-3n^l/2fj~l 

z=p^cP{3 + 2io)-y^H-^ 
wp = -KRle-^^+^P+H-^ 



= (l,-V3,0,0,0,±i/2) 
= (l,V3,0,0,±i/2,0) 

— VO|T:pg,=n 2(^:^+8) ' 2(^^+8)'^'^ 2(F+8) > 

= (1,^/3, 0,0,±V2,0) 

— V Oj^qig, zn 2(^2_,_8) , 2(€^+8) ' ^' ^ 2(F+8) > 
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The variables Wp and Wm are not independent one of the other and near an 
isotropic state, we have Wp oc Wm oc e~'^^H~^. Let us note moreover that 
W- is a positive variable. Three variables Wi are thus necessary to describe 
the curvature in the same way that three variables Ni are used in |^] for 
these Bianchi models. The Hamiltonian constraint is written as: 

xl + x'i + 2V + 12^2 + 12[w^{l + wi) ± 2w^{wmWpf/'^{l + w'i)+ 



w'^{Wm - 2wp)]{w'iw. 



and the field equations are: 

x+ = 72y'^x+ + 24{w^{x+ - 1)(1 + wt) ±W-{1 + 2x+){wmWpf/^{l + wl) 
+wl [{2 + x+)w'i-2{x+ -l)wl\}{w^^Wp)-^ (69) 
X- = 72i/2x_ + 2A{wl wtix- - \/3) + a;_ + Vs) ± W-{wmWpf/'^[w'i 

{-V3 + 2x-) + {y/3 + 2x_)] + w^_x- {wl^ - 2w^)}{w^_Wp)-^ (70) 
y = y{Uz + 72y'^ - 3 + 24[w':^{l + wi) ± 2{wmWp)^/'^W-{l + w'i)+ 



w_ 



{wl,-2wl)]{w'iwp) 1} 



i = y^(72z — 12£) + 24z[i/;p(l + wi) ± 2{wmWp 
W-iWm - 2wp)]{w'iwp)~^ 
Wp = Wp{-2 + 2x+ + 72y2 + 2A[w^{l + wi) ± 2w-{wmWp)^/^{l + u^? 
+w'i{wf^ - 2w^)]{w'iwp)^^} 



(71) 



(72) 



(73) 



w„ 



w, 



,{-2 - 2x+ + 72y2 + 24[wl{l + wi) ± 2w-{wmWp)^/^il + w^^) 



2 ^,,3 2wl)]{w'iwp)-^} 



+wt(w: 



W- 



2V3- 



W-X- 



(74) 
(75) 



The Hamiltonian equation for H becomes 
H -- 



H[72y^ + 24(±2^^i^ ± 2wy'wU'w. 



^K + ^) + i(7-2)A;2] 



(76) 



The signs it respectively represent the Bianchi type VIII or IX models. 
The constraint shows that the variables are not necessarily bounded: if one 
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of them diverges, this divergence can be compensated by that of Wm or Wp. 
Thus if we show that the isotropy occurs only for bounded values of Wm and 
Wp, it will mean that it occurs only for bounded values of all the variables 
(but (/)!). 

In order to achieve this goal, we will write that near an isotropic state Wp oc 
ii'„i w and W- — ^ const 7^ and to simplify the explanations which 
will follow, we will suppose that Wp Wm w and 1. Then the 

constraint of the Bianchi type VIII model shows that all the variables are 
positive and thus must take bounded values. With regard to the Bianchi 
type IX model, let us suppose that w diverges. Then if one puts x± = 0, 
we deduce from the constraint that 'iup' —> + 2:^ — 1/12 and from the 
equation for w that 'iuP' — > 3y^ — 1/12, implying that asymptotically — > 
and diverge like vo^ . However, with these limits one obtains from equations 
for y and z that y — ^ Mz^ — and i — > —12.^2;^ + 2z. Then the equilibrium 
for y and z can be only obtained when 2: ^ which is in contradiction with 
the divergence of z that we have just shown. One then deduces that an 
isotropic equilibrium state is impossible if Wp and Wm diverge and that all 
the variables must be bounded. It follows for the same reasons as for the 
previous Bianchi models, that the isotropisation is impossible for a bounded 
value of ri. 

One can also show that Wp and Wm can not tend to some non vanishing 
constants. Let us suppose that it is the case and define the two constants w 
and a such Wp ^ w and Wm aw. We introduce these limits into the 
equations for x± with x± = 0. It comes: 

x+ = -24u;2(l + it;_a3/2(l + wl) - 2w^_{l + a^) + wi)wZ'^ (77) 
x_ = -2AV3w^{w^ - 1)(1 - a3/2u;_ + w^)wZ^ (78) 

Then, for the Bianchi type VIII model, one deduces that equilibrium for x± 
will be reached only if q tends to the complex value (— 1)^/^ or/and if 
is negative what is impossible. For the Bianchi type IX model, equilibrium 
for x± can be reached if Wp Wm (i-e. /3± — > 0) and W- 1. Then, cal- 
culating the corresponding equilibrium points, the only ones which are real 
and such as Wp and Wm are different from are (x+, x^,y, z, Wp, Wm, W-) = 
(0, 0, ±{6i)-\{6i)-\±{l-e'^y/^{6£)-\ 1). They check the constraint equa- 
tion and are real if < 1. Moreover, one calculates that Wp and Wm tend 

to ±(1 - £2)^/2(1 - e e'^ + 36£^)~^/^ and thus reach equilibrium when 
Q — ^ -|-oo. Introducing these expressions into x+, it comes then that x+ 
tends to a complex value when ^ — -|-oo and thus that these equilibrium 
points are excluded. 

Consequently as for the previous models, the only possible isotropic equilib- 
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rium points are such as 



{x+,x^,y,z,Wp,Wm,W-.) = (0,0,±\/3 - £^{6V2) ^ ,i/6,0,0,W-o) 

w^o being a constant. The variables Wm and Wp behave asymptotically like 
^{i-e^)n gj^^ g2(i-£2)n follows that the asymptotic behaviors of the 

metric functions and the potential are the same ones asymptotically as for 
the other models. However, the sign of the Hamiltonian l(76|l is not preserved 
throughout the temporal evolution and it is thus not possible to know if the 
limit O — > — oo corresponds to the late or early times. 

4.5 Discussion 

Technically, compared to the Bianchi type / model, there exist several dif- 
ferences: 

• Put aside for the Bianchi type II and VIq models, the constraint does 
not imply directly that the variables x, y, z and w are necessarily 
bounded. It must be shown that it is the case when we approach a 
stable isotropic state. 

• It must be shown that the isotropy corresponds to a forever expansion 
of the Universe (17 — > —oo). 

• It must be shown that the product p±e^^ tends to zero. 

Physically, the models with curvature are more interesting than the models 
with flat space sections because they make it possible to show that then the 
isotropisation of class 1 is accompanied by an accelerated expansion and a 
flatness of the space sections. This comes owing to the fact that the equilib- 
rium points are such as the variables w related to the curvature disappear 
with the approach from isotropy, reducing the range of values in which the 
function i must asymptotically tend in order to allow the isotropisation. The 
asymptotic behaviors of the metric functions and the potential are then the 
same ones as for the Bianchi type / model because the Hamiltonian and the 
lapse function of the models with curvature behave asymptotically in the 
same way as for the flat model. Thus, we get the following result: 

Let us consider a minimally coupled and massive scalar-tensor theory and 
the quantity i defined by £ = ^(3^^,^)1/2 ■ asymptotic behavior of the 

scalar field to the approach of the isotropy is given by the asymptotic form of 
the solution for the differential equation cj) = 2 ^^_|^^^^ when 17 —00. This 
limit does not necessarily correspond to the late times for the Bianchi type 
VIII and IX models contrary to the other ones. A necessary condition for 
a class 1 isotropisation is that £^ < 1. If £ tends to a non vanishing constant, 
the metric functions tend to and the potential disappears like t~'^ . If £ 
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tends to vanish, the Universe tends to a De Sitter model and the potential to 
a constant. Whatever the curved Bianchi models, the Universe becomes flat 
and the expansion is accelerated. 

Thus, the accelerated behavior of the Universe and its flatness could find 
a natural explanation through the fact that it isotropises. Let us notice that 
the asymptotic behavior of the Bianchi type IX model is not oscillatory. 
However, this is not incompatible with a mixmaster behavior in the vicinity 
of a singularity as noted in PHI- The existence of a single isotropic equilib- 
rium state such as the curvature tends to vanish can appear shocking but it 
could be due to the fact that we apply the variability assumption of i. 
Considering again the scalar-tensor theory defined by — = \/2 and 

U = e™"^ that we studied in the previous sections, one thus finds that a 
necessary condition for isotropy will be m? < 2. A numerical integration for 
the case m = 0.8 and the Bianchi type IX model is illustrated on the figure 

m 



Figure 6: Evolution of the variables x±, y, z, Wp, Wm and «)_ when ^^'^^'^^ = 

U = e""*, i?Q = (\/247r)~^ and m = 0.8 with the initial values {x^, x- , y, z,Wp,Wm,w- , (/>) = 
(0.05,0.83,0.025,-0.12,0.02,0.002,0.2). All the variables tend to their theoretically predicted 
equilibrium values. 



5 Quintessence 

In this section, we will show that in the vicinity of an isotropic state, the 
scalar field may be quintessent under certain conditions. A quintessent scalar 
field is such that the ratio of its pressure Q and density (@J tends to a 
constant, this one having to be negative in order to be likely to explain the 
recent acceleration of the Universe. The scalar field pressure is then negative. 
To show the quintessence, one rewrites /j<^ and Ptf, using the variables of the 
Hamiltonian formalism. It comes: 

I^'/' + U/2 (79) 



2887r2i?6 
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In what follows, we will calculate for each Bianchi models the ratio w^p = 
P(t>/P4> and will determine the conditions such that the scalar field be quintessent. 

5.1 Bianchi type / model without a perfect fluid 

Using our previous results, we find that: 

= 2^ (81) 



and 



Then , we deduce: 



72t,'^RIxI 



2£ 3 2/^0 ,—2 
e oc t 



Thus with the approach of the isotropy, the scalar field behaves like a perfect 
fluid with an equation of state = w^p^ with = — 1 g [—1,1] and 
the function of the scalar field i can asymptotically be interpreted like the 
barotropic index of this fluid. It will be quintessent if i'^ < 3/2 what is 
compatible with the isotropisation. The WMAP^ data show that < 
—0.78 what gives us < o.33 and a deceleration parameter q = — 1 < 
-0.67 



5.2 Bianchi models with curvature without a perfect fluid 

The scalar field energy density and pressure behave in the same way but the 
barotropic index is this time in the range w^f) G [-1,-1/3]: the scalar field 
is always asymptotically quintessent when the Universe isotropises. It is in 
agreement with the fact that its expansion is always accelerated. 



5.3 Bianchi type / model with a perfect fluid such as 

The only thing which changes when we approach an isotropic state compared 
to the section WTT] is that < Calculations relating to the scalar field 
energy density and pressure are thus the same ones as in the absence of a 
perfect fiuid but from now on tt;^ S [— 1,7 — 1]: according to the value of 7 
the scalar field will be (if 7 < 1) or not (if 7 > 1) always quintessent. 

'http://lambda.gsfc.nasa.gov/ 
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5.4 Bianchi type / model with a perfect fluid such as fi^ -h 

In this case we have U oi pm what results naturally in: 

U = - (X pm (82) 

If the scalar field is such as p^f, oc p^, that means that its state equation must 
be the same one as that of the perfect fluid because it is the only means of 
explaining why their energy densities, which are obtained by writing their 
conservation laws, scale in the same way. Let us check it. Using the limits 
above, we calculate that: 

and 

^ 7(2-7) 3^0 

Consequently, from H79l8n|l . we obtain that P41 = {"f — ^)pcf> , which is well in 

agreement with (|H^ and the fact that to the approach of the isotropy, the 

2 

metric functions tend to t^'' as if there were only the perfect fluid and not 
the scalar field: obviously, the scalar field is never quintessent. 

5.5 Conclusion 

When we consider a scalar field, we found that the constant to which tends 
the function £ at the time of the isotropisation can asymptotically be in- 
terpreted as being the barotropic index of the state equation characterizing 
the scalar field. This one can then be quintessent with the approach of the 
isotropy for the Bianchi type / model with or without a perfect fluid such as 
— > if £2 < 3/2. Without surprise, quintessence is impossible if 7^ 0. 
On the other hand, in presence of curvature it is always the case because 
isotropisation imposes that i'^ < 1. 

6 General conclusion and perspectives 

In this chapter, we considered the properties of the homogeneous cosmolog- 
ical models in scalar-tensor theories with a minimally coupled and massive 
scalar field and studied their isotropisation process. 

From the point of view of dynamical systems analysis, we detected three 
families of equilibrium points, corresponding to three different ways for the 
Universe to reach isotropy and called them class 1, 2 and 3. We were inter- 
ested in class 1 and obtained results consisting of: 

1. the location of the stable isotropic equilibrium points 
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2. some necessary conditions for their existence, constraining the scalar- 
tensor theories 

3. the asymptotic behaviors of the metric functions and scalar field. 

Those of these results related to asymptotical behaviors were calculated un- 
der the assumption that the Universe tends sufficiently quickly to its equi- 
librium state. Mathematically speaking that means that when we approach 
the isotropy, the variables appearing in the field equations tend sufficiently 
quickly to their equilibrium values so that one can neglect their variations 
in the calculations. We showed how this assumption could be raised with 
regard to the function i. On the other hand for the other variables, a per- 
turbative study to the approach of equilibrium could be necessary and some 
progress will have to be made in this direction to supplement the study of 
the class 1 isotropisation process. 

Finally, the state in which the Universe is when it reaches the isotropy show 
some interesting characteristics which can be summarized in the following 
way: 

• the function i = (j)U^U~^{3 + 2ijj)~^l'^ must tend to a constant whose 
biggest value is constrained by the presence of curvature or perfect 
fluid. 

• the Universe is expanding such as the metric functions tend to power or 
exponential proper time laws and the potential respectively disappears 
like or tends to a constant. 

• the presence of curvature favor late time acceleration and quintessence. 

• the Universe is asymptotically flat. 

For the Bianchi type / model, our results can be summarized in table |21 In 
presence of curvature but without a perfect fluid, the results are similar but t 
must be smaller than one, thus supporting the appearance of an accelerated 
expansion and of quintessence. 

Other work is in hand on this subject. Thus in jSHI; for the Bianchi type / 
model when the scalar field is not minimally coupled, i.e. when we consider 
a gravitation constant varying with the scalar field, we have also manage to 
constrain the scalar-tensor theories in such a way that they are compatible 
with the isotropisation. The asymptotic behaviors of the metric functions 
and the potential were also given but in the Einstein frame where 4> is min- 
imally coupled. It is indeed impossible to obtain these behaviors without 
using quadrature when the field is not minimally coupled. A case of class 2 
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Isotropisation 


Asymptotical 


Quintessence 




necessary conditions 


behaviors 






e^<3 


if f ^ 0, e-^ ^ t^"' 

and U 
if ^2 ^ 0, e-^ ^ e* 
and U — > const 


= (If - l)p^: 
Yes iff < 3/2 




£2 < f 


idem as above 


idem as above 






and U 


P</. = (7 - 1)^0: 
No 



Table 2: isotropisation of the Bianchi type / model with a minimally coupled and massive 
scalar field when the Universe reach sufficiently quickly its isotropic state. 



isotropisation was detected. 

In fSnj, always for the Bianchi type I model, we generalized the results ex- 
posed in this chapter in the presence of two scalar fields. Such theories can 
be related to compactification process or to the presence of complex scalar 
fields. In this last case, we have numerically observed that isotropisation was 
mainly of class 3. 

In [HO], we considered the Bianchi models with curvature in the presence of a 
perfect fluid such as 7 G [1,2] and showed that, contrary to the Bianchi type 
/ model, the results obtained in the absence of perfect fluid were unchanged 
when Qrn whereas isotropisation was impossible when Qrn A except 
if 7 < 2/3. 

Finally in 001, we showed an interesting link between the scalar-tensor the- 
ories leading to isotropisation and those allowing a flattening of the rotation 
curves for the spiral galaxies: in both cases, the function i must tend to a 
constant and the scalar field can be quintessent. 

We hope having defined an operational framework able to guide future re- 
search on the subject of the isotropisation of homogeneous but anisotropic 
cosmologies in scalar-tensor theories. Amongst other things, it would be 
significant to take into account the cases of a negative potential or Brans- 
Dicke coupling function. It would then correspond to a non respect of the 
weak energy condition, assumption which often appears in the literature but 
which still misses physical motivations. More significant, it would be use- 
ful to generalize our results when the convergence to the isotropic state is 
not sufficiently fast to neglect the terms of the second order for the integral 
of (variability assumption of i) or for the bounded variables {xi,y, z,Wj) 
which were used to rewrite the Hamiltonian equations. Lastly, an important 
extension of our results would consist of the study of classes 2 and 3 isotropi- 
sation which were approached only numerically through some applications 
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in papers that we evoke above. 
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